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Abstract 

By means of density functional calculations the magnetic and electronic properties and phase stabilities of the Heusler compounds Co2MSi 
(with M=Ti,V,Cr,Mn,Fe,Co,Ni) were investigated. Based on the calculated results we predict the ferromagnetic phases of the compounds 
Co2TiSi, Co2VSi and Co2CrSi to be half-metals. Of particular interest is Co2CrSi because of its high density of majority spin states at Fermi 
energy in combination with a reasonably high estimated Curie temperature of 747K. The compounds Co2TiSi and Co2VSi are thermodynam- 
ically stable, whereas Co2CrSi is a metastable phase which might be stabilized by suitable experimental techniques. 



I. INTRODUCTION 



n. COMPUTATIONAL DETAILS 



In the pioneering work of de Groot et. al.' NiMnSb 
and PtMnSb were predicted to be half-metal ferromagnets, 
for which the majority-spin states are of metallic character 
whereas the minority-spin states have a gap at Fermi en- 
ergy. A larger class of Heusler alloys has such peculiar elec- 
tronic and magnetic properties^, which -in combination with 
large magnetic moments and high Curie temperatures- makes 
these materials attractive for the design of single-spin electron 
sources^ and spin injectors^ in the field of magnetoelectron- 
ics and related technological applications. 

Within the series of Co2MSi compounds (for 
M=Ti,V,Cr,Mn,Fe,Co,Ni), CoiFeSi and CoaMnSi were 
studied to some extent. For Co2MnSi a Curie temperature of 
Tc = 985 K and a total magnetic moment of about 5 fiB^'^'^ 
were measured. For Co2FeSi, the measured Curie temper- 
ature of Tc = 1100 K is the highest for all known Heusler 
alloys reflecting the large total magnetic moment of about 6 
HB^. Recently, thin films of CozMnSi and CoaMnGsii^ 
were fabricated, for which their magnetic properties are still 
under debate. Density functional theory (DFT) calculations 
for Co2MnSi^ '^ confirmed its half-metal properties and the 
measured large magnetic moment. A conventional DFT study 
on Co2FeSi''^ derived that Fermi energy does not fall any 
more into the gap of the minority-spin states. Treating the 
d-like states of Fe as strongly localized states, the half-metal 
ferromagnetic property could be recaptured, as it is claimed 
by experiment»iS 

In the present work, by means of a DFT approach 
we examined the series of Heusler alloys Co2MSi (with 
M=Ti, V,Cr,Mn,Fe,Co,Ni) assuming they crystallize in the typ- 
ical L2i structure. For all these compounds we derived mag- 
netic and electronic properties as well as thermodynamical 
stabilities and elastic properties. In the seven compounds, 
only for M=Ti, V, Mn, Fe compounds with the L2i struc- 
ture were synthesized. ^-^-'^.iSjie ^j^g experimentally claimed 
structure of Co2CoSi (or CosSi) is of the DOig type at tem- 
peratures of about 1200 °C whereas below 1192 °C the com- 
pound could not be stabilized.'^ Concerning Co2CrSi and 
Co2NiSi, no experimental data are available in literature. 
Nevertheless, for the sake of completeness and comparison, 
Co2CrSi, Co2CoSi, and Co2NiSi with the L2i structure are 
included in our study. 



Our DFT calculations were performed by application of the 
plane wave Vienna Ab initio Simulation Package (VASP) — 
in its projector augmented wave formulation for the poten- 
tials. The exchange-correlation potential and energy were de- 
scribed within the generalized gradient approximation (GGA) 
of Perdew and Wang'^ in combination with the parametriza- 
tion of Vosko et al.'" for spin polarized densities. The cu- 
bic lattice parameters were optimized self-consistently, and 
care was taken to converge the total energy in terms of basis 
functions and of k points for the Brillouin zone integration. 
Ferromagnetic as well as some selected antiferromagnetic or- 
derings were considered. Local properties (such as projected 
densities of states and local magnetic moments) were deter- 
mined within suitably chosen spheres centered at the atomic 
positions. Finally, thermodynamic and elastic stabilities were 
derived by a procedure reported in Ref. |2ll For Co2CrSi 
and Co2FeSi, additional studies were made by treating the d- 
like states as strongly localized states in terms of a so-called 
LDAh-U approach^^ with only the one effective parameter U- 
J, similar to the LDAh-U approach of Ref. 10. Our LDAh-U 
calculations were done on the basis of the GGA potentials as 
used for the conventional calculations. 



III. RESULTS AND DISCUSSIONS 

FiglU summarizes trends of the density of states (DOS) 
and related quantities which were derived from the de- 
scribed VASP calculations applying GGA potentials. For 
M=Ti,V,Cr,Mn the Fermi energy Ep lies in the pronounced 
gap of the minority-spin states (left side panels), determining 
the half-metal character of these compounds. According to 
FigHJ' the width of the band gap has values in the range of 
0.6 to 0.7 eV . The Heusler alloy Co2FeSi, however, is not a 
half-metal anymore, because Ep falls into an uprising peak of 
the minority-spin states. The -now indirect- band gap is only 
0. 1 eV wide, being strongly reduced due to a flat tail of the 
mentioned peak (see corresponding left side panel). For the 
half-metal compounds, however, this pronounced peak is to- 
tally unoccupied. For M=Co,Ni the gap vanishes, because all 
the d-bands are nearly filled and the separation of the Co-d and 
M-d band center is now too small for the formation of a gap. 
The crucial gap for the half-metal compounds, as discussed by 
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FIG. 1: (Color online) Magnetic and electronic properties derived 
from VASP calculations with GGA potentials for the compounds 
Co2MSi (M=Ti,V,Cr,Mn,Fe,Co,Ni) with L2i crystal structure. Pan- 
els on the left side: ferromagnetic spin-polarized density of states 
split into majority-spin and minority-spin states. Panels on the right 
side: panel a) shows N(Ep), the density of states at the Fermi level 
for majority-spin states (full line) and minority-spin states (broken 
line); panel b) shows Eg (in eV), the characteristic band gap at Ep of 
the minority-spin states; panel c) illustrates the trend of m (in fin), 
the local magnetic moments of Co (broken line) and M atoms (full 
line), together with experimental data (crosses). 



Galanakis et al. for Co2MnSi^, is due to a strong hybridiza- 
tion between Co-d and M-d states, combined with large local 
magnetic moments and a sizeable separation of the d-like band 
centers. We found that a related strong hybridization feature 
(a small gap for M=Ti,V or a deep valley for M=Cr,Mn,Fe) 
occurs already in the non-spin polarized DOS, for which Ep 
cuts through a pronounced peak indicating instabiliy. When 
allowing for ferromagnetic spin polarization the Fermi energy 
is then pinned in the gap of the minority-spin DOS as long 
as there are not too many states to be filled and the gap is 
sufficiently large, which is the case for the compounds with 
M=Ti,V,Cr,Mn. In these cases, the local magnetic moments 
m]y[ of the M-atoms increase linearly, whereas ruQ^ remains 
rather constant (see FigH};). It should be noted, that the nu- 
merical values of the local moments depend on the choice 
of the localization criterion (e.g. on the choice of radii of 
atomic spheres in our case), and therefore -in general- they 
are not integer numbers. However, the total magnetic mo- 
ments rritot are exactly integer numbers for the true half-metal 



TABLE I: Results of VASP calculations with GGA potentials: 
magnetic moments (total moment mtot in /^s/f.u., local moment 
niQ^^ and mjyj in /is), density of states at Fermi level (N(|,Ep) 
for majority spin and N(J,,Ep) for minority spin states, in states 
eV~^ spin~^ atom~^ ), band gaps Eg (in eV) of the minority 
spin state, and estimated Curie temperatures Tc (in K) for Co2MSi 
(M=Ti,V,Cr,Mn,Fe,Co,Ni). 



Co2TiSi Co2VSi Co2CrSi Co2MnSi Co2FeSi 



"^Co 


1.01 


1.03 


0.98 


1.02 


1.34 


myi 


-0.05 


0.82 


2.08 


2.99 


2.79 


miot 


2 


3 


4 


5 


5.48 


N(T,-Bf) 


0.24 


0.35 


0.80 


0.31 


0.18 
















0.85 




0.62 


0.68 


0.72 


0.66 


0.10 


Tc 


385 


566 


747 


928 


1109 



TABLE II: Calculated bulk moduli and elastic constants (cn, C12, 
C44 in GPa) for the CoaMSi (M= Ti, V, Cr, Mn, Fe). 



M 


Ti 


V 


Cr 


Mn 


Fe 


B 


215 


216 


227 


221 


204 


Cll 


303 


255 


297 


316 


247 


C12 


172 


197 


193 


174 


182 


C44 


126 


130 


145 


143 


133 


c' 


65 


29 


52 


71 


33 



compounds, as shown by Table U 

Comparing Co2MnSi to Co2FeSi, the local moment mpg 
is now smaller than m^^^, and the linear trend for m-^ (for 
M=Ti,V,Cr,Mn) is now destroyed, as shown by Fig[T]; . The 
reason is that according to a conventional DFT GGA calcu- 
lation Co2FeSi is not a half metal. This is consistent with 
a recent DFT study by Wurmehl et. alM^, who also found 
that the Fermi level is not in the band gap of the minority- 
spin DOS when applying GGA or the local density approxi- 
mation (LDA) for the exchange-correlation potential. Treat- 
ing the d-states as strongly localized states by means of an 
LDAh-U approach and choosing suitable parameters for U-J^, 
the compound Co2FeSi became a half-metal with a local mo- 
ment of 6 HB in agreement to experiment. We reproduced 
the LDAh-U results which can also be obtained from conven- 
tional GGA calculations for a volume which is expanded by 
5% compared to the calculated equilibrium value. It should 
also be mentioned that in the case of Co2MnSi, our results 
agree well with earlier DFT studiesi^ii^ 

According to our results, amongst the three predicted half 
metals Co2MSi (M=Ti,V,Cr) the Cr compound is the most in- 
teresting one, because of its large density of states of the ma- 
jority spins at Fermi energy of N(t,Ep) = 0.80 states eV^^ 
spin^^ atom^^. This is the largest value for all known fer- 
romagnetic half-metals. The reason for this large value is 
that the Fermi energy cuts through strongly localized states 
of mostly Cr-d like character, as illustrated by the band struc- 
ture and d-like DOS in Fig.|2] The contribution of Co d-states 
to N(t,Ep) is very small, because Ep falls into a deep min- 
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TABLE III: Calculated lattice parameter (a in A ) and enthalpies of 
formation (AH, kJ (mol of atoms)^^) compared to available experi- 
mental data for Co2MSi (M= Ti, V, Cr, Mn, Fe). NM: non magnetic 
calculations, FM: ferromagnetic calculations. 



M 


Ti 


V 


Cr 


Mn 


Fe 


> 


a NM 


5.7205 


5.6393 


5.5897 


5.5582 


5.5410 




FM 


5.7609 


5.6621 


5.6295 


5.6457 


5.6231 


Li. 

LiJ 


exp. 


5.770° 


5.659 * 




5.654 


5.640 


LiJ 


AH NM 


-61.5 


-32.2 


-11.7 


-9.42 


-4.7 




FM 


-64.4 


-40.7 


-29.7 


-44.9 


-33.9 





"Ref[Tl 
''Ref [l4lfT^ 



imum of as well as states. According to this finding, 
Co2CrSi is a half-metal which offers 100% spin-selectivity 
combined with a pronounced site selectivity of strongly local- 
ized Cr states at E^?. This is in contrast to the well-known 
half-metals Co2MnSi and Co2FeSi, for which the states at Ef 
are of more delocalized character. 

The Curie temperatures Tc for Co2MSi (M=Ti,V,Cr) were 
estimated similar to a model presented in Ref. J^by applying 
the relation Tc = 23 + ISlmtot, which is linear in the total 
magnetic moment rritot per unit cell. The results are 385, 566, 
and 747 K for Co2TiSi, CozVSi, and CojCrSi, respectively. 
Again, the Cr compound is of special interest, now because 
of its rather large T^, making it interesting for technological 
applications. 

Because the Cr d-states at E^ are strongly localized one 
might argue that these states have to be treated in a more suit- 
able way e.g. in terms of the LDAh-U approach as described 
previously for the study of Co2FeSi. In our case, we applied 
the approach of Dudarev et al. — which needs only the dif- 
ference A = U ~ J of the on-site Coulomb and exchange 
parameters. Similar to the study Wurmehl et a/.'" we chose 
A ~ 4.8 eV for Co, whereas for Cr A was varied between 
and 3 eV. (The choice of A = corresponds to the conven- 
tional GGA calculation.) For A = 2 eV, Co2CrSi is still a 
half-metal with a strongly increased gap for the minority spin 
states of 1.9 eV However, for A = 3 eV the Fermi energy 
cuts through minority spin states, and the half-metal charac- 
ter is destroyed. The main effects of the localization enforced 
by the LDAh-U treatment are an increased gap width and the 
occupation of more majority spin states than for the conven- 
tional GGA or LDA approaches which lowers the Fermi en- 
ergy. This lowering effect could be so strong (i.e. for larger 
A), that finally E p cuts through minority spin states below the 
gap. The applied LDAh-U approach has the advantage, that 
only one unknown parameter needs to be introduced, namely 
the difference of U minus J. Most of the other LDAh-U ap- 
proaches (e.g. Refl23l) rely on the two independent parameters 
U and J. A suitable choice of these parameters is -however- at 
present not possible, because no experimental information is 
yet available. For Co2CrSi being a half-metal or not, more 
elaborate calculations might be necessary to determine either 








C\ m 


n. 





































































W K W L 



W K 



(a) band structure 




-4-2 2 

E-Ep (eV) 



6 -4 -2 
E-Ep (eV) 



(Z?) projected DOS 



FIG. 2: (Color online) VASP results with GGA potentials for ferro- 
magnetic CoaCrSi: spin polarized band structure (upper panel) and 
projected density of states (lower panel). 



the unknown parameters, or to go beyond the limitations of 
the LDA+U approaches. 

Measurements of the magnetic and electronic properties of 
Co2CrSi would give the answer. However, the preparation of 
suitable samples is not straightforward, because Co2CrSi is 
thermodynamically metastable. (In contrast to Co2TiSi and 
Co2VSi, which are even known to exist thermodynamically.) 
According to our ab initio calculations of the enthalpy of for- 
mation, the compound Co2CrSi is unstable against a thermo- 
dynamical separation into the phases Co2Si, Cr and Si by 
3 kJ (mol of atoms)^^. On the other hand, Co2CrSi is sta- 
ble against decompositions into CoaSi, Cr and Si (by 23 kJ 
(mol of atoms)-!), q^^^^ Qj. ^^^^ gj^j^^ ^j^^j Qj. 

(by 29.7 kJ). Obviously, Co2CrSi might be stabilized by ad- 
vanced experimental techniques for synthesizing metastable 
states. Elastically, however, the Cr compound is very sta- 
ble as demonstrated by the calculated elastic constants of 
cii=297 GPa, Ci2=193 GPa, and C44=145 GPa, and the bulk 
modulus of B=227 GPa (see Table|llli. 
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The calculated lattice parameters listed in Table |III] for 
the known Ti-, V-, Mn-, and Fe-based half-metals are -as 
expected- in good agreement to experiment 7'"'''"*i'^''^ In this 
table, enthalpies of formations are also given. In particular, 
non spinpolarized (NM) calcualtions are compared to ferro- 
magnetic (FM) calculations, which result in significant energy 
gains, in particular for the Cr-, Fe-, and Mn-based alloys. The 
calculated elastic constants are given by Table HH which show 
no anomolous behaviour In particular, the Co2CrSi alloy fits 
in the trend of the other compounds, and is certainly elasti- 
cally stable. (For Co2CoSi and Co2NiSi with their assumed 
artificial structures the derived shear moduli c' are negative, 
indicating elastic instability). 



Co2VSi, and Co2CrSi. The Cr compound is of particular in- 
terest because of its high density of states of 100% polarized 
states at Fermi energy together with a strong selectivity of Cr- 
d states and an elevated Curie temperature. We hope that our 
work stimulates experiments on Co2CrSi once this metastable 
phase could be synthesized. 



IV. CONCLUSION 

Summarizing, based on density functional calculations 
we predict three ferromagnetic half-metals, namely Co2TiSi, 
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